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Abstract 

In a method for approximate solution of Sturm-Liouville equations and related spectral problems 
was presented based on the construction of the Delsarte transmutation operators. The problem of 
numerical approximation of solutions and eigendata was reduced to approximation of a primitive of 
the potential by a finite linear combination of certain specially constructed functions obtained from the 
generalized wave polynomials introduced in [15], [26]. The method allows one to compute both lower 
and higher eigendata with an extreme accuracy. 

Since the solution of the approximation problem is the main step in the application of the method, 
the properties of the system of functions involved are of primary interest. In |27| two basic properties 
were established: the completeness in appropriate functional spaces and the linear independence. In this 
paper we present a considerably more complete study of the systems of functions. We establish their 
relation with another linear differential second-order equation, find out certain operations (in a sense, 
generalized derivatives and antiderivatives) which allow us to generate the next such function from a 
previous one. We obtain the uniqueness of the coefficients of expansions in terms of such functions and a 
corresponding generalized Taylor theorem, as well as formulas for exact expansion coefficients involving 
the operations mentioned above. We also construct the invertible integral operators transforming powers 
of the independent variable into the functions under consideration and establish their commutation 
relations with differential operators. We present some error bounds for the solution of the approximation 
problem depending on the smoothness of the potential and show that these error bounds are close to 
optimal in order. Also, we provide a rigorous justification of the alternative formulation of the proposed 
method allowing one to make use of the known initial values of the solutions at the left endpoint of the 
spectral problem. 


1 Introduction 

One of the important mathematical tools for studying problems related to Sturm-Liouville equations was 
introduced in 1938 by J. Delsarte [S] and called nni the transmutation operator. It relates two linear 
differential operators and allows one to transform a more complicated equation into a simpler one. Nowadays 
the transmutation operator is widely used in the theory of linear differential equations (see, e.g., DP: 01: m, 
[30], [3T], [33]). Very often in the literature the transmutation operators are called the transformation 
operators. It is well known that under certain regularity conditions the transmutation operator transmuting 
the operator A = + q{x) into B = can be realized in the form of a Volterra integral operator 

with good properties. Its integral kernel can be obtained as a solution of a certain Goursat problem for 
the Klein-Gordon equation with a variable coefficient. In spite of their attractive properties and importance 
there exist very few examples of the transmutation kernels available in a closed form (see [23]). 

Several recent results concerning the transmutation operators made it possible in m, [19] to convert 
the transmutation operators from a purely theoretical tool into an efficient practical method for solving 
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Sturm-Liouville equations and related spectral problems. The method is called the analytic approximation 
of transmutation operators. It is based on the result from m where a new complete system of solutions for 
the Klein-Gordon equation was constructed. The functions of that system are called generalized wave poly¬ 
nomials. Since the integral kernel of the transmutation operator is a solution of that Klein-Gordon equation, 
it is possible to approximate the integral kernel as well. It was shown in m how the problem of approxima¬ 
tion of the integral kernel reduces to the solution of two one-dimensional problems of approximation of the 
pair of functions giix) = ^ + j fg q{s)ds and g 2 {x) = j lo <l{s)ds in terms of specially constructed families 
of functions {c„} and {s„}, appearing as traces of the generalized wave polynomials on the line x = t. With 
respect to the variable of integration the approximate kernel results to be a polynomial. This is especially 
convenient since to obtain solutions of the Schrodinger equation the transmutation operator is applied to the 
functions sin ^/Xt and cos V^t, solutions of the simplest such equation Bv = Xv, and thus all the involved 
integrals can be calculated explicitly. One of the advantages of the method is that due to the independence 
of the integral kernel of the spectral parameter the error of the computed eigendata does not increase for 
higher eigenvalues. One can compute, e.g., the 1000th eigenvalue and eigenfunction with roughly the same 
accuracy as the first ones. 

In [27| we presented a rigorous justification of the method, however several important questions remained 
unanswered. It was observed in various numerical examples that the method allows one to obtain highly 
accurate eigendata and demonstrates exponential convergence with respect to the number of the functions 
{c„} and {sn,} used while we only proved that the analytic approximations converge without any convergence 
rate estimates. Also the functionality of the method was proved for the symmetric segment [— 6 , b] and in 
such form the method did not allow one to make use of the known initial values of the solutions at a; = 0 . 

The present paper is dedicated to further study of the analytic approximations of transmutation oper¬ 
ators and of the systems of functions {c„} and {s„} involved. In particular, we show that the systems of 
functions {c„} and {s„} are closely related with another linear differential second-order equation, find out 
certain operations (in a sense, generalized derivatives) which allow us to present Taylor-type formulas in 
terms of the linear combinations of the functions {c„} and {s„}. We also construct the invertible integral 
operators transforming powers of the independent variable into the functions {c„} and {s„} and establish 
their commutation relations with differential operators. We prove that the convergence rate of the analytic 
approximations method depends on the smoothness of the potential and provide the justification of the 
alternative formulation of the method allowing one to make use of the known initial values of the solutions 
at an endpoint of the spectral problem. 

The paper is structured as follows. In Section [5] we introduce some necessary notations, definitions and 
properties concerning special systems of functions called formal powers, generalized wave polynomials and 
their traces {c„} and {s„}. We present the definition of the transmutation operators, recall some properties 
and briefly outline the method of analytic approximation of the transmutation operators. Also we present 
new relations for the integral kernel of the transmutation operator on the characteristics x = t and x = —t. 
We would like to mention that these relations (12.181) and ()2.20p already found applications to a posteriori 
accuracy verification in In Section [3] we show that the functions {c„} and {s„} are closely related to 
the formal powers for the equation y" — q{x)y = 2Ay'. In Section |4] we introduce the generalized derivatives 
7 i and 72 acting as 7271 c„ = ns„_i and 7271 s„ = nc„_i (the second formula is valid for n>2 only). Also 
we study Taylor-type formulas in terms of the functions {c„} and {s„}. We present the formulas for the 
coefficients and prove the corresponding Taylor-type theorem with the Peano form of the remainder term. 
In Section [5] we obtain a new representation for the preimage of the integral kernel of the transmutation 
operator. Also we prove that the even coefficients of the generalized Taylor series for the two functions related 
to the integral kernel coincide. In Section | 6 ] we study the Goursat-to-Goursat transmutation operators Gi 
and G 2 mapping the powers of the independent variable to the functions c„ and s„ respectively. We prove 
some commutation relations involving these operators, derivatives and generalized derivatives 7271 - We 
show that the existence of certain number of the generalized derivatives ( 7271 )^ of a function implies a 
certain order differentiability of the preimage under the action of either operator Gi or operator G 2 and 
vice versa. In Section [7] we show that there exists a close relation between the smoothness of the potential q 
and the convergence rate of the analytic approximation of the transmutation operator, corresponding direct 
and inverse theorems are proved. In Section |S] we present alternative proofs of the main theorems of the 
analytic approximation method. The proofs are based on the well-posedness of the Goursat problem and 
does not involve neither inverse operators nor pseudoanalytic function theory. Also the proofs justify the 
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applicability of the analytic approximation method in the restricted setting, when the equation in considered 
on the half-segment [0, b] only and no continuation of the potential onto the whole segment [—&, b] is used. 


2 Transmutation operators and systems of functions {c„} and {s„} 


2.1 Formal powers 

Let / € C[a, b] be a complex valued function and f{x) ^ 0 for any x € [o, b]. The interval (a, b) is assumed 
being finite. Consider two sequences of recursive integrals 

= 1, X^'^\x)=n f ds, Xq € [a,b], n = l,2,... 

-'Xq 

and 


= X^^'>{x)=nf (/2(s))^ ' ds, Xo€[a,b], n = l,2,.... 

J Xq 

Define two families of functions {^kYk^Q and constructed according to the rules 

t j k odd, 

^f{x)X^^\x), k even. 


( 2 . 1 ) 


and 


'ipkix) 


xW(a;) 

fix) 

X^’^\x) 

fix) 


k odd, 
k even. 


( 2 . 2 ) 


2.2 SPPS representations 

The following result obtained in m (for additional details and simpler proof see [18] and [22]) establishes 
the relation of the system of functions and Sturm-Liouville equation. 

Theorem 2.1. Let q be a continuous complex valued function of an independent real variable x S [a, 1>] and 
X be an arbitrary complex number. Let f be a solution of the equation 

f" -qf = 0 (2.3) 

on (a, b) such that f G C^ia, b) fl b] and f{x) 0 for any x G [a, b] (it always exists, see E.em.ark \2.S\) . 
Then the general solution y G C^ia, b) fl C^[a, b] of the equation 

y” -qy = Xy (2.4) 

on (a, b) has the form y = Ciyi + C 22/2 where Ci and C 2 are arbitrary complex constants, 

00 ^k °° 

and = (2.5) 

and both series converge uniformly on [a, b] together with the series of the first derivatives which have the 
form 



f 

yV^2fe + 2k'4’2k 



and ^2 = E 




fe=0 


(2A: + 1)! 


f 

—jT^k+l + (2/c -f 1) ■i/'2fc 


The series of the second derivatives converge uniformly on any segment [ai,6i] C ia,b). 


( 2 . 6 ) 
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Representations (1^ and (12.61) are called SPPS (Spectral Parameter Power Series) representations. They 
have been implemented for solving a variety of spectral and scattering problems related to Sturm-Liouville 
equations. The first work using Theorem 12.11 for numerical solution was [22] and later on the SPPS method 
was used in a number of publications (see 0,0,0, m, [S], m, BS], m, bhi and references therein). 

Remark 2.2. It is easy to see that hy definition the solutions yi and 2/2 from (12.51) satisfy the following 
initial conditions 


yiixo) = fixo), y'lixo) = f'{xo), 

y2ixo) = 0 , ^2(2^0) = l//(a;o)- 

Remark 2.3. It is worth mentioning that in the regular case the existence and construction of the required f 
presents no difficulty. Indeed, let q be real valued and continuous on [a, 6]. Then (12.31) possesses two linearly 
independent real-valued solutions fi and /2 whose zeros alternate. Thus, one may choose / = /i + 1/2 • 
Moreover, for the construction of f\ and /2 in fact the same SPPS method may be used In the case of 
complex-valued coefficients the existence of a non-vanishing solution was shown in \2S[ Remark 5], see also 

m- 


2.3 Transmutation operators 

Let if be a linear topological space and Ei its linear subspace (not necessarily closed). Let A and B be 
linear operators: Ei ^ E. 

Definition 2.4. A linear invertible operator T defined on the whole E such that Ei is invariant under the 
action of T is called a transmutation operator for the pair of operators A and B if it fulfills the following 
two conditions. 


f. Both the operator T and its inverse T ^ are continuous in E; 
2. The following operator equality is valid 


AT = TB 


(2.7) 


or which is the same 


A = TBT 


-1 


Our main interest concerns the situation when A = + qix), B = and q is a continuous 

complex-valued function. Consider the space E = C[—b, b] with b being a positive number. In 0 and [23] 
a parametrized family of transmutation operators for A and B was studied. Operators of this family can 
be realized in the form of a Volterra integral operator with a kernel associated to a value of the complex 
parameter h, 

pX 

( 2 . 8 ) 


Tu(x) = u{x)-\- j 'K{x,t\h)u{t)dt. 

J —X 


ft is convenient to interpret the parameter h as follows. Let / be a solution of (12.31) satisfying the initial 
conditions 

/(O) = 1 and /'(O) = h. 

Then there exists a unique operator of the form (12.81) satisfying (12.71) and such that T [1] = /. We will denote 
it by T/ and its kernel by 'K{x,t). Thus, the operator T/ has the form 


Tfu{x) = u{x)-\- f 'K{x,t)u{t)dt. 

J —X 


(2.9) 


The kernel can be dehned as K{x,t) = H(^^, |t| < \x\ < &, H being the unique solution of the 

Goursat problem 

^^^^^=q{u + v)li{u,v), ( 2 . 10 ) 


q{s)ds, 


H(0,u) = | 


du dv 

H(..,0) = i + 1 


0 
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( 2 . 11 ) 



where h := /'(O). 

If the potential q is continuously differentiable, the kernel K itself is a solution of the Goursat problem 

( 2 . 12 ) 

(2.13) 


/ 5 ^ \ 

K(a:,a:) =-^ + qis)ds, K{x,-x) = f^. 


If the potential g is n times continuously differentiable, the kernel li.{x,t) is n + 1 times continuously 
differentiable with respect to both independent variables. 

The following theorem states that the operators T f are indeed transmutations in the sense of Definition 


(2.14) 


Theorem 2.5 ([26]). Let q G C[—b,b]. Then the operator Tf defined by (12.911 satisfies the equality 

d^ 



d^ 

j T/M = Tf 



for any u G C'^[—b,b]. 

Remark 2.6. Tf maps a solution v of the equation v" + = 0, where to is a complex number, into a 

solution u of the equation 

u” — q{x)u + uj^u = 0 (2-15) 

with the following correspondence of the initial values u{0) = v(0), u'(0) = w'(0) + hv(0). 

The integral kernel K(a::,t) satisfies in the region |t| < |x| < & the following integral equation (see [261 
(3.3) and (3.4)] and [30l §2]), 


K(*,() = i + 5 


q{s) ds + 


q{a + /3)K(a + fi,a — /3) d/3 da. 


(2.16) 


Differentiating ()2.16|1 with respect to t we obtain, c.f., |2T1 (3.11)], that 

K 2 (a;, t) = ^q ^ q{z)K{z, x + t- z)dz-^ q{z)K{z, z - {x - t)) dz, (2.17) 

here and below by Ki(a;,t) and K 2 (a;,t) we denote the partial derivatives with respect to the first and the 
second variable respectively. 

Let us introduce the following useful notation 


Q{x) := / q{s)ds. 
Jo 


Proposition 2.7. The following relations for the partial derivatives of the transmutation kernel are valid 
on the line t = x, 

If,. . . Q‘^{x)' 


and 


Ki(a;, x) = - U(a;) + hQ{x) + ^ 


K 2 (a;,a;) = ^ ^q(a;) - hQfix) - 9.-^ 


(2.18) 

(2.19) 


Proof. We obtain from (12.171) and (I2.13P that 

K2{x,x) = p{x) q{z)K{z,z)dz=^q{x)-^J^ + q{s) ds'^ dz 


which coincides with (12.191) observing that q(z) q(s) ds dz = Q^(x)/2. 
Let us notice that 

Ki(x,x) + K 2 (x,x) = = ■^<?(a:). 


dx 


Hence (|2.18p immediately follows from (|2.19l) . 


□ 
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Similarly we obtain the following proposition. 

Proposition 2.8. The following relations for the partial derivatives of the transmutation kernel are valid 
on the line t = —x, 

Ki(a;, -x) = K 2 (a;, -x) = + ^Q{x). (2.20) 

Relations (I2.18P and (|2.20p were already utilized in to estimate the error of the approximation of the 
derivatives of solutions of equation (IQ) . 

As can be seen from (I2J]) . the definition of the transmutation operator requires the knowledge of its 
integral kernel K in the region |t| < |a;| < b. However the integral kernel K is well-defined (via (I2.10|) ~ (I2.11I) ') 
in the larger region S : |a:| < b, |t| < b and is continuously differentiable (twice for q G C^[—6, b]) there, see 
P5] . For the rest of this paper we consider the integral kernel K to be defined in this larger domain S, it 
allows us to present many results in a simpler and more natural form. First of such results is the following 
proposition defining the inverse operator 

Proposition 2.9 l[23jl. The inverse operator T^. ^ can be represented as the Volterra integral operator 

Tf^u{x) = u{x) — J 'K{t,x)u{t) dt. 


Together with the transmutation T/ it is often convenient to consider another couple of operators enjoying 
the transmutation property (12.141) on subclasses of C^[—b, b] (as well as on subclasses of 6]), for details 
see m and additionally |24], 

Tcw{x) = w{x) + / C{x,t)w{t)dt 
Jo 


and 


Tsw{x)=w(x)+ f S{x,t)w{t)dt 

Jo 


with the kernels C and S related to the kernel K by the equalities 


C{x, t) = K.{x, t) -1- K(a:, —t) 


and 

S(x, t) = K(a;, t) — K(x, —t). 

The following statement is valid. 

Theorem 2.10 ([30]). Solutions c(uj,x) and s{uj,x) of equation (12.151) satisfying the initial conditions 


c(cc,0) = l, c'Ju;,0) = h (2.21) 

s(a;,0) = 0, 4(w,0) = l (2.22) 


can be represented in the form 

r 

c{uj,x) = cosojx + / C{x,t)cosujtdt 
Jo 


and 


s{uj, x) = 


sinwcc 

OJ 


, sin wt , 
S(a;, t) - dt. 


(jj 


The following important mapping property is a corollary of Theorem 12.11 


Theorem 2.11 ([3], [33]1. Let q be a continuous complex valued function of an independent real variable 
X G [—b,b] and f be a particular solution of (12.31) such that f ^ 0 on [—b,b] and normalized as /(O) = 1. 
Denote h := f'{0) G C. Then 

Tf [a;^] = ipk{x) for any k G Nq. (2.23) 
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Remark 2.12. The mapping property (12.231) of the transmutation operator reveals that the SPPS represen¬ 
tations (1231) from Theorem \2.1\ are nothing but the images of Taylor expansions of the functions cosh \/Xx 
and sinh '‘/Xx under the action of T j. 

In what follows we assume that / 0 on [—b,b], /(O) = 1 and denote h := f'(0) G C. 

We introduce the following two systems of functions 


Cm (a;) = 

even k—0 
m 

Sm{x) = 


odd k—1 


m 


x^ipm-k{x), TO =1,2,... and Co{x) = uo{x,x) = f{x), (2.24) 


x^ipm-k{x), TO = 1,2,... and So = 0. (2.25) 


Example 2.13. In a special case when / = 1 we obtain that Cm{x) = Sm{x) = 2"* to = 1, 2,_ 

The systems of functions {cml^^o {®m}m=o 8'i'ose in [2^ as traces of so-called generalized wave poly¬ 
nomials which form a complete system of solutions of (12.121) and are used in [26] for uniform approximation 
of the transmutation kernels K, C and S. 

Definition 2.14 ([E]). The following functions 

™ /to\ ™ /to\ 

Uo=(poix), U2m-lix,t) = Y Y]T 7 n-k{x)t'‘, U2m{x,t)= Y Y j Tjn-k{x)t'", 

even k—0 ^ '' odd k—1 '' 

are called generalized wave polynomials (the wave polynomials are introduced below, in Example 12.151) . The 
following parity relations hold for the generalized wave polynomials. 

Uo(x,-t) = Uo{x,t), U 2 n-lix,-t) = U 2 „-lix,t), U 2 n{x,-t) = -U 2 n{x ,t) . 

Example 2.15. In a special case when / = 1 we obtain that (pk{x) = k G Nq and Uk{x, t) = pk{x, t) where 
Pk are wave polynomials [m Proposition 1] defined by the equalities 

Po{x,t) = l, P 2 m-lix,t) = Y f + 

even k—0 '' ' 

m . \ . 

{x,t)= Y + 

odd fe=l ^ ' 


P2 

Remark 2.16. From Theorem \2.11\ one obtains that = Ty [pn]- 

Recall that S denotes a closed square on the plane (x, t) with a diagonal joining the endpoints (6, b) and 
(-5,-6). 

Theorem 2.17 ([27]). Let the complex numbers ao, - ■ -,gat and hi,... ^bN be such that 


h 1 
2+4 


/ q{s)ds - Y. a-nCn{x) 

•^0 n=0 


< £l 


and 


\ C ^ 

- q{s)ds - Y, Ksnjx) 
a 0 „_1 


< ^2 


for every x € [—b,b]. Then the kernel K(x,t) is approximated by the function 


N 


N 


Kn{x, t) = aoUo(x, t) + Y anU 2 n-l{x, t) +Y ^nU 2 n{x, t) 

n—1 n—1 

in such a way that for every (x, t) G S the inequality holds 

|K(x,t) — KN{x,t)\ < e 

where e > 0 depends on ei, £2 and q. 


(2.26) 


(2.27) 


(2.28) 
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Remark 2.18 ([H]). If the functions ^ ^ q{s)ds and j q{s)ds admit uniformly convergent on [—6, b] 

respective series expansions 

^ nx OO ^ OO 

2 4 / ^ X! 4 / d(s)ds = ^ bnSn{x), 

® n=0 n—1 


then the kernel K admits a uniformly convergent in S series representation 

OO OO 

K{x,t) = aoUo{x,t) + ^ anU2n-l{x,t) + ^ bnU2n{x,t). (2.29) 

n=l n=l 

In [27) the uniform approximation of the kernel from Theorem 12 .1 71 was applied to obtain a result on the 
uniform approximation (with respect to x and ui) of the solutions c(oj,x) and s(a;,a;) of equation (12.1511 . 

Theorem 2.19 ([27]). The solutions c{uj,x) and s{uj,x) of equation p.l5p satisfying (12.211) and (|2.22|) 
respectively can be approximated by the functions 


^ ^ fn\ 

C]s[{uj, x) = coswa: + 2 H ( \ipn-k{x) I t^ COS Ujtdt 

n—0 even k—0 ^ ^ ^ 


(2.30) 


and 


SAr(a;, x) = — ( sinwa: + 2 { ]ipn-k{x) [ t^sinujtdt] 

odtt=l ) 


(2.31) 


where the coefficients and {&»}„=! axe the same as in Theorem |il.77| and the following estimates 

hold 

, . , , .. £sinh(C'x) 

|c(w,x) - CAr(w,x)| < 


and 


|s(w,x) - SAr(a;,x)| < 


C 

£sinh(Cx) 

\oj\C 


for any w £ C, w 7 ^ 0 belonging to the strip |Ima;| < C, C > 0, where e > 0 depends on ei, £2 and q. 

The existence of the appropriate number N and coefficients {a„}^^Q and {bn}n=i required in (12.261) and 
(12.271) is established by the following statement, see also Section [7| 


Proposition 2.20 ([27]). The systems of functions {c„})yQ and {s„})yj^ are linearly independent and 
complete in C^[—b,b] and Cg[— 6 , &], respectively, with respect to the maximum norm. Here Cq [— 6 , 6 ] denotes 
a subspace of C^[—b,b] consisting of functions vanishing in the origin. 

In the next section we show that there exists another way for construction of the functions {c„})(yQ and 
{s„}^i which additionally reveals some of their properties. 


3 Functions {c„} and {s„} as formal powers 

We shall use the following result from [^. Consider the Sturm-Liouville equation of the form 

N 

(p(x)it')'+ q(x)it = [u] , X e {a,b) (3.1) 

where Rk are linear differential operators of the first order, Rk [u] := rk(x)u + Sk(x)u', k = 1 , ... N, and the 
complex-valued functions p, q, Vk, Sk are continuous on the finite segment [a, b]. 
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Define the formal powers for equation as follows 


^{-n) ^ ^ Q forn G N, 

_^( 0 ) ^ ^( 0 ) ^ 




X^^\x) 


N 


f{s)J2Rk /(s)A’("-2'=+l)(s) 
ds 


k=l 


ds, 


^ ^Xn 


>xo Pis)p{s) 

' Xq 

ox. N 


P (s)p(s)’ 

ps)Y,Rk [/(s)A’("-2'=+i)(s)j ds, 


fc=l 


n - odd, 


n - even, 


n - odd. 


n - even 


where / is a particular complex-valued solution of equation m for A = 0 and xq is an arbitrary point of 
the segment [a, b] such that p{xo) ^ 0. 

Theorem 3.1 (SPPS representations for polynomial pencils of operators, [IS]). Assume that on a finite 
interval [a,b], the equation 

{p{x)v'y + q{x)v = 0 

possesses a particular solution f such that the functions fRk[f], k = and are continuous on 

[a, b]. Then the general solution of (13.11) has the form u = ciUi +C 2 U 2 , where ci and C 2 are arbitrary complex 
constants and 

00 00 

ui=fY^ and U 2 =fY^ yn;t"( 2 n+i) ^ (3 2 ) 

n—0 n—0 

Both series in (|3.2I1 converge uniformly on [o, b]. 

Below it becomes clear that of particular interest is the equation 


y" - q{x)y = 2Aj/' 


(3.3) 


where g is a continuous complex-valued function on [0,6]. Its general solution can be obtained according to 
Theorem 13.11 

Corollary 3.2. Assume that the function f is a solution of the equation 

v" — q(x)v = 0 

on (0,6) such that /(O) = 1 and f{x) 7^ 0, a; G [0,6] (see B.emark 1 2. S\ according to the existence of such a 
nonvanishing solution). Then a general solution of (13.31) can be written in the form 


yi=fJ2 and y2 = fJ2 

n—0 n—0 

with the formal powers and defined as follows 

y(0) ^ y(0) ^ 

2/o'/(«) ds, n - odd, 


Y^^\x) = 




ds 


n - even. 


Pis)' 

'2/o'/(s) (/(s)T(”"^Hs))'rfs> n - even, 
= { . ds ,, 

1 fo ^ (s) j^> n - odd. 


(3.4) 


9 










This is a direct corollary of Theorem 13.II where N = 1, Ri \y\ = 2y', etc. 

From now on we assume that / satisfies the conditions of Corollarv l3.2l 

Remark 3.3. It is easy to verify that the solutions yi and y 2 from (13.41) fulfil the following initial conditions 

2/i(o) = i, y;(o) = /'(o), 

2 / 2 ( 0 ) = 0 , 2 / 2 ( 0 ) = !. 


In the following proposition a relation between the functions c„ and s„ defined in (12.241) and (12.251) with 
the formal powers from Corollarv l3.2l is established. 

Proposition 3.4. 


= jy(2n-i)^ ^ = / ^ gdd n G N 

nl n\ \ / 

and 

n! n! V / 

Proof. Observe that m is a solution of the equation 

u” — q{x)u = }fu 


(3.5) 

(3.6) 


(3.7) 


iff z; = is a solution of 

v" — q{x)v = 2\v'. 

Indeed, for v = we have v' = Xv + e^^u' and 

v" = Xv' + Xe^'''u' + = Xv' -\- X{v' — Xv) + [qu + }?u) = 2Xv' + qv. 

Analogously, the function w = e~^^u solves the equation w" — q{x)w = —2Xw'. 

Now let us consider the solution u of (1X71) having the form u = ui + U 2 with the solutions ui and U 2 
constructed according to Theorem 12.II 




fc=0 


We have then 
Consequently, 


f(0) = 1 and u{0) = f'{0) + A. 


(3.8) 


V = 


OO ■. h, U OO A 

A'^x'^ X^ 


fe=o ■ i=o 


g>k-j 


j\ (k — ?)! 

fc=0 j=0 ' k=0 


OO 

i:^‘()t+S)- 


In a similar way we obtain that 


w = e = 


OO 

EA^(t-S)- 


(3.9) 


(3.10) 


fe=0 


On the other hand, the functions v and w can be constructed according to Corollary 13.21 For this we 
observe that 


^(0)=u(0) = l, z;'(0) =/'(O) + 2A, 

w;(0) = 1, n;'(0) = /'(O). 


Comparing these values with those from Remark 13.31 we obtain 


V = yi+ 2Xy2 


f [ A”r(2") + 2 ^ a”+if(2 


n+l) 




n=0 


(3.11) 


10 










and 


(3.12) 


zz; = /^(-A)”F( 2 n)^ 

n—0 

Comparing (13.91) with (13.111) and (13.101) with (I3.12p we arrive at the relations 

oo oo 

E t + S) = ^ ^^ 

k—0 k—1 

and 

OO OO 

k=0 ■ ■ k=0 

Adding and substracting these two equalities, due to the uniform convergence of the power series with respect 
to A, we obtain (13.51) and (13.6p . □ 


4 Generalized derivatives 


Following the ideas from m. [m, m. let us introduce the following couple of operations which will be 
called the generalized derivatives, 


7i5(a;) 


fi^)- 




dx \f{x) 


and 


i2g{x) 


1 r g'{s)ds 

2 Jo /(s) 


The following proposition clarifies their relation to equation (1531) . 


Proposition 4.1. For any g G C^[0,6], 


727i5(a;) = \ g{s)ds. 

Moreover, the operator 7271 can be extended onto [0, b] by the rule 

72715 ( 2 :) = J (^g'{x) - g'{0) - q{s)g{s) ds^. 

Indeed, by definition, 


(4.1) 


( 7 ®))= U’ 


/o fis)ds\ ds\f{s)JJ 2 Jo f{s)ds 
^ ^ ■{f{s)g''{s) - f'{s)g{s))ds = J j {g"{s) - q{s)g{s))ds 


{fis)g'{s) - f{s)g{s))ds 


lo /(s) 


= ^(g'ix)-g'{0)- q{s)g{s)ds^. 
Thus, for a solution of (13.31) we have 


Note that 


72712 /( 2 :) = X{y{x) - 1 /( 0 )). 


71 / = 0 


and 


7271 [1] (a:) 


1 

2 



q{s)ds 
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Remark 4.2. Below we consider the functions ( 7271 )^ [1], j = 0,1,.... They appear in IttOl Lemma 1.4-1] 
where it was shown that they are well defined for j < n if q G and the function ( 7271 )" [ 1 ] 

possesses a square integrahle derivative. Moreover, it can he easily verified using (14.11) that if q G 6 ] 

then the functions ( 7271 )^ [ 1 ] are well-defined (and are continuous functions) for j < n + 1 . 

In what follows when considering the function ( 7271 )^ [1], if not specified explicitly, we suppose that 
qGWi-^[0,b]. 

In the following statement we summarize several properties of the functions c„ and s„ related to the 
generalized derivatives. 


Proposition 4.3. The following relations are valid 


Proof. Consider 


7 iCo = 0, 71C1 = 1, 71S1 = /^, 

7271C0 = 7271C1 = 0, 7271S1 = / - 1, 

7iCn(0) = 7is„(0) = 0, n = 2,3,..., 
7271 ^^ — ^ — 1; 2, ... , 


7271 = nc„-i, 


("-!)! 

Cn^i 


( 7271 )'’ -7 = ■^ 

n! 


= ("-!)! 

(7271^” [/ - 1] 

0 


("-!)! 


( 7271 )'^ ^ j) 

n\ 


(7271^” [/ - 1] 

0 


n = 2,3,..., 

if j is odd and j < n, 
if j is even and j < n, 
if j P n > 0 and n is even, 
otherwise, 

if j is odd and j < n, 
if j is even and j < n, 
if j P n and n is odd, 
otherwise. 


7isi = f y(i)y ^ y(l) I 

Here we used (13.51) and the definition of the formal powers and . Note that 


(4.2) 

(4.3) 

(4.4) 

(4.5) 

(4.6) 


(4.7) 


(4.8) 


yW(a;)=2/ f{s)fis)ds = fix)-l. 
Jo 


Hence 71 S 1 = /^. The other two equalities in (14.21) can are proved in a similar way. 
Equalities (14.41) follow from Proposition 13.41 and the definition of and Y^^\ 
Equalities (14.31) follow trivially from (14.21) . 

For n = 2, 3,... consider 


Cn 1 

7271 —T = o727l 

n\ 2 


/(?■<“"' 


This is according to Proposition [TH Hence 

f ((1 + (-1)”) + 2y(^"-i)y 

(1 + (-1)") 


Cyj 1 

7271—j- = 7:72 
n! 2 

1 

= 272 


1+ (-!)' 


(^f{s)Y^^^-^\s))' ds+ f (/(s)y(2»-3)(5)yrf5 


fix) (^i7_LiZ!y(2"-2)(2.) +y(2n-3)(2.)^ 
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fix) ((1 - (-1)”-^) + 2 r( 2 (»-i)-i)(a;)^ 


1 
2 

Sn—1 
(n- 1)!' 


Equality (14.611 can be proved in a similar way. 

The proof of (j4.7l) and (14.81) follows from the preceding equalities. 


□ 


Definition 4.4. Functions of the form 


N 


Cn = Y1 


0!.r 


n—0 

where an are complex numbers, will be called c-polynomials of order and functions of the form 

N 


Sn = 


n—1 


n t 5 

n\ 


where /3„ are complex numbers, will be called s-polynomials of order N. 

We stress that in general Cn and S^r are not, of course, polynomials. 

Proposition 4.5. Let Cn be a c-polynomial of order N then its coefficients are uniquely determined by the 
relations ao = CAr(O) and 


3 

^i+i = 7i (7271^ CAr(O) + y] an7i ( 7271 )'^”" [1] (0), j = 0,1, 2,..., - 1. (4.9) 

even n—2 

Proof. Assume j < N being odd. Then due to (j4.7l) . 

N , . 

il 2 liy Cn{x) = ^an (7271 )■' 

= £ a„(727i)""”[/-l](a^)+ Y. 

even n—2 n—j-\-l ^ 

Now applying to this equality 71 and considering it at cc = 0 we obtain 

3 N 

7i ( 7271 )^ Cat(0) = Y «n7i ( 7271 )'^”" [/- 1] (0) + Y “’"■ 7 ^ 3 'br 

even n—2 
3 

= Y «ri 7 i (7271)'^”" [/- 1 ] ( 0 ) + aj+i- 

even n—2 

Finally, observing that 71 / = 0 we obtain (|4.9I) . The proof of (14.91) for any even j is similar. The recursive 
relation (|4.9p implies the uniqueness of the coefficients q;„. □ 

Analogously the statement concerning s-polynomials is proved. 

Proposition 4.6. Let S^r be an s-polynomial of order N then its coefficients are uniquely determined by 
the relations 

3 

/ 3 i-n = 7i (7271 )'’Sat( 0)-f Y ^n 7 i (7271)'^”" [1] (0), j = 0,1,2,..., V - 1. 

odd n—1 
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Lemma 4.7. Let the complex valued function F of the variable x be such that in the neighborhood of x = 0 
its generalized derivatives 71 ( 7271 )'^ F exist for j = 0,... N — 1 and are continuous at x = 0. Let be a 
c-polynomial with the coefficients of the form ao = ^( 0 ) and 

j 

a^'+i = 7i ( 7271 )'^-P(O) + X! «n7i ( 7271 )’’”" [1] (0), j = 0,1,2,.... (4.10) 

even n—2 

Then for the function R := F — C^v the following relations are valid 

R{Q) = 7 ii?( 0 ) = ... = 71 ( 727 i)'^”' ^(0) = 0. (4.11) 

Proof. Obviously, i?(0) = F(0) —CAr(O) = F(0) —ao = 0 and 7 ii?( 0 ) = 7 ii^( 0 ) — 7 iCAr( 0 ) = 7 iF( 0 ) —ai = 0. 
For any j = 1,... — 1 we obtain 

7i ( 7271 )^ -R(O) = 7i ( 7271 )'’ FiO) - 7i ( 7271 )'’ Cn{0) 

3 

= 71 ( 7271 )'’-P'(O) - aj+i + ^ an 7 i ( 7271 )'^”” [ 1 ] ( 0 ) = 0 

even n—2 

due to (14.101) . □ 

A similar statement is true for s-polynomials. 

Lemma 4.8. Let the complex valued function F of the variable x be such that in the neighborhood of x = 0 
there exist 71 ( 7271 )'^ F for j = 0,... N —1 and are continuous at x = 0, F(0) =0. Let Sjv be an s-polynomial 
with the coefficients of the form 

3 

/3,+i =7i(727i)'>(0)+ ^ /3„7 i ( 7271 )'”” [1] (0), j = 0,1, 2,..., iV - 1. (4.12) 

odd n—1 

Then for the function R := F — S^v the following relations are valid 

R{0) = 7ii?(0) = ... = 71 (727i)'^”' ^(0) = 0. 

Lemma 4.9. Suppose the complex valued function R of the variable x be such that in the neighborhood of 
X = 0 there exist 71 ( 7271 )'^ R for j = 0,... iV — 1 and are continuous at x = 0. Let 

i?(0) = 7 ii?( 0 ) = ... = 71 ( 7271 )” -R(O) = 0. 

Then R{x) = o(a:"'+^), a: —>■ 0. 

Proof. The proof is by induction. For n = 0 by assumption we have that i?(0) = 7 ii?( 0 ) = 0. Note that 
by definition of 71 , jiR{0) = —f'{0)R{0) + R'{0) (we took into account that /(O) = 1), and hence from the 
assumption we obtain that R'{0) = 0 which implies that R(x) = o{x). 

Now, assuming that the statement is true for n—1, let us prove it for n. For this, consider Ri := 7271 !?. 
Then by assumption we have 7 ii?i( 0 ) = ... = 71 ( 7271 )"”^ i?i(0) = 0 as well as i?i(0) = 0 due to the 
definition of 72 . Hence Ri{x) = o{x^). Also we observe that 

Ri{x) = \ - q{s)R{s)) rfs = ^ - J q{s)R{s)ds'^ , (4.13) 

where we took into account that R'{0) = 0. 

Due to the mean-value theorem we have 

R{x) = R{x) — R{0) = X (Rei?'(ci) -I- i Imii'(c 2 )), 

where Ci and C 2 are some points between 0 and x. Using 1)4.13p we obtain 

2i?i(ci)-|- J q{s)R{s)ds'^ + ixlm ^Ri{c 2 ) + J q{s)R{s)ds 

= X (Re (o(x"') -I- o{x'^~^^)) -I -1 Im (o(a;") -I- o(x"“'"^))) = o{x'^~^^). 

□ 


R{x) = xRe 
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Theorem 4.10 (Taylor-type theorem with the Peano form of the remainder term). Let the complex valued 
function F of the variable x be such that in the neighborhood of x = 0 its generalized derivatives 71 ( 7271 Y F 
exist for j = 0,... N — 1 and are continuous at x = 0. Then 

N . . 

F{x) = + o{x^) (4.14) 

n\ 

n—O 

where the coefficients a„ are defined as in Lemma [f71\ 

If additionally F(0) = 0, then also 

nl 

n—l 

where the coefficients fin are defined as in Lemma lff^ 

Proof. For the proof of (14.141) one needs to observe that due to Lemma HTTl R := F — '^n=o satisfies 
(14.lip and hence according to Lemma |49l R{x) = o{x^), a; —>■ 0. 

Equality (14.1511 is proved similarly with the aid of Lemmas 14.81 and 14.91 □ 

Remark 4.11. Consider the solution u of (IT7P introduced in the proof of Proposition \3.4\ as satisfying 
dSH). Then the function cosh Xx ■ u{x) is an example of a function admitting a uniformly convergent on [0, b] 
Taylor-type expansion in terms of the functions c„. Indeed, we have 

00 / \ 

cosh Xx ■ u[x) = ^ A" ^ 


and analogously, 

00 / \ 

sinh Ax • u{x) = A” ^" . 

^ ^ n\ 

n—l 

The relations are obtained by adding and subtracting the functions v and w from the proof of Proposition 

[M 

Example 4.12. Lemmas 14.71 and 14.81 easily give ns the exact expressions for several first coefficients a„ and 
fin of the generalized Taylor formulas for the functions Fi{x) = f + 7 Jq q{s) ds and F 2 (x) = j q{s) ds 
appearing in (12.261) . (12.271) and in the representation (12.291) . see also (|6.241) . First, observe that for a function 
g such that g{0) = 0 one has 7 ig( 0 ) = g'(0). Second, we can compute several expressions ( 27271 )^( 1 ] using 
(HID- We have 


( 2727 i)[ 1 ](x) = -Q{x), 

( 2727 i)^[ 1 ](x) = -q{x) + ^(O) - ^Q^{x), 

(27271 )^[ 1 ](a;) = -q'ix) -f q'{0) - {q{x) -f g( 0 ))Q(x) -f J q'^{s) ds -f ^Q^{x), 

(27271 )‘^[l](x) = -q'fx) + q"{0) - q'{x)Q{x) - q{0)q{x) 

+ ^^(O) + ^q(x)Q^(x) - g(s)( 2727 i)^[l](s) ds. 

Note that by the definition of the operator 7271 all the expressions above are equal to zero at x = 0. Hence, 

7 i[l]( 0 ) = -h, 7 i( 2727 i)[ 1 ]( 0 ) =-g( 0 ), 71 ( 27271 )^ [1] ( 0 ) =-g'( 0 ), 

71 (27271 )"[1](0) = -g"(0), 71 (27271 )"[1](0) = -g"'(0) - 2q'i0)q{0). 
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For Lemma K7\ the function -F’i(x) = ^ — ^( 727 i)[l], hence we obtain 


9(0) 


ao = 2 ’ '^1 = - 2 ^ *^2 = 


9 '(0) hq{Q) 


03 — 


9"(0) _ V(0) fe^9(0) 

16 4 4 ’ 


04 = 


4 ’ 

q"'{Q) h{q"{0)-2q\Q)) 


32 


16 


For Lemma IDHl the function F 2 {x) = — 5(7271 )[1], hence we obtain 


/3i = 


m 

4 ’ 


P2 = 


9 '(0) hq{0) 


o 9 "( 0 ) 9 ^( 0 ) 

* = ^6 - 8 “' * = 


4 ’ 

g"'(0) /i(g"(0)-292(0)) 


32 
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Observe that the coefficients with the even indices coincide. Later in Propositions 15.31 and 16.91 we give the 
proof of this phenomenon. 


5 The preimage of the transmutation kernel 

Together with the function K let us consider its preimage k := ^ [K]. Due to Theorem [23] it is a solution 

of the wave equation and hence admits the representation 


k(a;,r) = 


- V' 


where the functions (p and ip are unique up to an additive constant. 
Proposition 5.1. The function k(a;,r) admits the representation 


k(a;,r) = ^ + i^ q{s)ds - J dt 


K t 


K U,-- 


dt. 


Proof. From Proposition 12.91 we have 

k(x, t) = K(a;, r) — f 'K{t,x)K.{t,T)dt. 

J —X 

From this equality and (EH) we obtain 

fx h I P" 

h{x,x) = p {x) + ip{0) = K{x,x) — / K'^{t,x)dt = — + - / q{s)ds — / K'^{t,x)dt 

J-x 2 2 Jq J 


(5.1) 


(5.2) 


where (12.131) was used. 
Similarly, 


/■“ h r 

\i{x,—x) = (p{0)+ip{x)='K{x,—x)— / 'K{t,x)'K{t,—x)dt — — — / K.{t,x)'K{t,—x)dt. 

J—X ^ J—X 


Note that 


Thus, 




k(0,0)= 9^(0)+1/^(0) = -. 




(5.3) 

(5.4) 
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and 



Adding these expressions and taking into account (15.41) we obtain (15.21) . 
Remark 5.2. The function ip' is even. Indeed, due to (lOl) . 


dt — ip (0). 


□ 


ip{—x) = — — p{Q)— / K.{t,—x)'K.{t,x)dt 

^ Jx 

h 

= — — (/? (0) + / K(t, x)'K{t, —x)dt = —ip{x) + h — 2p{Q). 

^ J-X 

From this fact and from it follows that 

ki(a:, r) = ki(r, a:) and k 2 (a;, r) = k 2 (r, a;). 

Proposition 5.3. Suppose that the conditions from Remark \2.18\ are fulfilled and hence K admits the 
representation (12.291) . Then 

a 2 n = b 2 n for any n S N. (5-5) 

Proof. Indeed, in this case 

OO 

k(a;, r) = aoPo{x, t) + X! {o.nP2n-i{x, t) + bnP2n{x, r)) 

n—1 

h 1 ^ 

= 2 + 2 X! ■ 

n—1 

Thus, 

^ OO 

= (0) + 2 X! 

n—1 

and 

^ (“4“) = ^ ^ 

From the last equality we have 

1 “ 

Ip {x) = Ip {Q) F 2 ” (a„ - 6 „) x" 

n—1 

and 

^ OO 

if (-x) = -(/; ( 0 ) + - ^ (- 1 )” 2 " (a„ - bn) x”. 

n—1 

Due to Remark 15.21 we obtain the equality 



^(-ir 2 " (a„- 6 „)x" = 

n—1 

which proves the equality of the even coefficients (1531) . 


^2" (a„-6„)x" 

n—1 


□ 
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6 Goursat-to-Goursat transmutation operators and the general¬ 
ized derivatives 


By S we denote a closed square with a diagonal joining the endpoints {b,b) and {—b,—b) (c.f., Sect. 

Let O := d'^ — df and the functions u and u be solutions of the equations Dm = 0 and (□ — ( 7 ( 1 )) m = 0 in 
S, respectively such that u = TfU. In [37] the operator G acting on the space C^[—b,b] x Co[— 6 , 6 ] as 


^ 1 / u{x, x) + u{x, —x) 

2 \u{x,x)—u{x,—x) 


1 /'u{x,x) + u{x,—x)\ 

2 \u{x, x) — u{x, —x) J 


( 6 . 1 ) 


was introduced. It is bounded together with its inverse and is related to the operator transforming the 
Goursat data for the equation Dm = 0 into the Goursat data for the equation (□ — (/(a;))?/ = 0. The 
operator G allowed us to prove the completeness of the functions {c„})(^q in G^[—b, b] and of the functions 
{s„}^i in Cg)— 6 , 6 ], see [271 Corollary 4.5]. The following proposition summarizes some properties of the 
operator G. 

Proposition 6.1 ([37]). 

1) The operator G defined by (]6.1I) on C^[—b,b] x Cg]—6, 6 ] admits the following representation 


G 


T]{x) 


Gi[rj{x)] 

G2K(a^)] 


G, 


r]{x) - 

G.iax) 


+ 


v(0) 


where G+ and G- have the form 


G±r]{x) = T]{x) + J K(a;, t) (^r] ± V 


t — X 


dt. 


2) Both operators G+ and G_ preserve the value of the function in the origin and Gj. : G^\—b,b] —>■ 
C^[-b,b], G_ :Ggi[-6,6]^Ggi[-6,5]. 

3) There exist the inverse operators Gf^ and GZ^ as well as the inverse operators Gf^ and Gf^ defined 
on C^[—b,b], and the inverse operator for G admits the representation 


G“ 


-1 ( vix ) 

^(x) 




V{x) - ^ 

GZ^ [fix) 


+ 


r)(0)' 


4) The operator G maps the powers of x into the functions c„ and s„ according to the following relations 


G : 


1 -^ 


Co (a;) 
0 


G : 2 


1-1 X 


0 




Cn(a^) 

0 


G: 2’^ 


0 




0 

s„(a;) 


Together with the operators G+ and G_ consider the following pair of operators acting on C^[—b,b] 

't — x' 


and 


r_? 7 (a;) = ?7(0) + f K{x,t)r] 
J —X 

T+r]{x) = f]{x) + f K{x,t)f] 

J —X 


t + X 


dt 


dt. 


Proposition 6.2. The images of the powers of x under the action of the operators r± are given by the 
following relations 

r_[a;”](a:) = ^(c„(x) - s„(x)), r+[a;”](a;) = ^(c„(a;) + s„(a;)), n > 0. (6.2) 
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Proof. We present the proof for the first relation in (16.21) . the proof for the second relation is similar. 
If n = 0, then 

r_[l](a;) = l+ f K(a;,t) • Idt = T/[l](a;) = Co(a:) = Co(a;) - So(a;). 

J — X 

If n > 0, then using Remark l2.16l and formulas from Example 12.151 we obtain 


T_[x^]{x)= j K{x,t) 
J — X 


t-x\ 1 r 

— -P2nit,x)) dt 


1 1 r 

= ^{P 2 n-l{x,x) -P 2 n{x,x)) + — J K{x, t) {p 2 n-l{t, x) - P 2 n{t, x)) dt 
= ■^{u 2 n-iix,x) - U 2 n{x,x)) = ^(c„(a;) - s„(a;)). 


□ 


Proposition 6.3. The following commutation relations hold for the operators r_ and r_|_ on C'^[—b, b], 


{d^ - q{x))r+= dT+d and {d^ - q{x))T_ =-dT_d, (6.3) 

moreover, on C^[—b,b] the operators r± satisfy 

2727ir+= (r_|_ — i5)i9 and 2^2'yi^-=— d)d, (6.4) 

where d denotes d/dx and S[r]\(x) = r]{0). 

Proof. Note that the relations (j6.3l) and (16.41) are equivalent on C^[—b,b] due to Proposition 14.11 

Relations (lOD can be verified directly similarly to [231 Theorem 6]. We omit this lengthy verification 
and exploit Proposition 16.21 to verify (lOl) . 

Let n > 1. Then using Propositions 14.31 and 16.21 


2727ir±[a:"] = 27271 


1 

2 n 


(Cn 


■^^{sn ± c„) = ±r±[nai” = ±(r± - (5)[9a;”], 


i.e., (16.41) holds for the functions a:" for any n > 1. The validity of (16.41) for the cases n = 0,1 can be easily 
verified as well using Proposition 14.31 By linearity (16.41) holds for any polynomial. 

Consider a function p{x) G C^[—b,b] and let Pn be a sequence of polynomials such that Pn{0) = ??(0), 
Pn( 6 ) = v'i^) s-nd p'ni^) d'i.x) uniformly as n ^ 00 . Then necessarily Pn(x) —>■ ri{x) uniformly for 

X G [—b, b]. Using Proposition 14.11 the relations (16.41) for the polynomials can be written as 

-^{^+[pr^{x)]) =p'„(o) + /ip„(o) + ^ ( 7 (s)r+[p„(s)]ds + r+[p;(x)]-p;(o) 

and ^ 

^(r_[p„(a;)]) = hpniO) + J g(s)r_[p„(s)] ds - r_[p(j(a;)] +p(i( 0 ) 

The right-hand sides of these expressions converge uniformly as n —> 00 . The differentiation operators in the 
left-hand sides are closed operators (defined on suitable subspaces of C^[—b, b]) hence allowing us to pass to 
the limit as n —>■ 00 and verify (16.41) for the function 77 . □ 

The operators Gi and G 2 have the following representations in terms of the operators Pi and r 2 , 

Gi=T++T_ - f-S, and G 2 = r+ - P- -b <5. (6.5) 

Applying the generalized derivatives 7271 to (16.5p and using (16.41) we obtain the following corollary. 

Corollary 6.4. The following commutation relations hold for the operators Gi and G 2 on G^[—b,b], 

27271 Gi = (G 2 -(5)5 and 27271 (G 2 - 5) = (Gi+ (/-2)5)5. ( 6 . 6 ) 
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Repeated application of Corollarv l6.4l allows us to obtain the following result. 

Corollary 6.5. Suppose that q G C'^[— 6 , b] and let a function u G C"‘[— 6 , b] be given. Then there exist all the 
functions ( 7271 )-^ [Giu] for j < min{m,p + 4} and all the functions ( 727 i)^ [G 2 M — u(0)] for j < min{TO,p + 3} 
and they are given by the following expressions. 

If n = 2k then 


n 

(2727irGiu = Giu(")+/-u(")(0)-2 ^ u«(0) • (27271 )”■'[!], (6.7) 

even j—2 
n—1 

(2727i)"[G2M-u(0)] = G2u(”)-u(")(0)-2 ^ . (27271 (6.8) 

odd j—1 

If n = 2k + 1 then 

n— 1 

(2727i)”GiM = G 2 u(”)- u(”H0)-2 ^ (27271 )"-J'[l], (6.9) 

even j—2 
n 

(2727i)”[G2M-u(0)] = Giu(")+/. J”)(0)-2 ^ u(^)( 0)(27271 )”-^[l]. (6.10) 

odd j—1 

Proof. Applying (16.6p twice we obtain the following commutation relations valid on G^[—5, 6 ], 

( 2727 i)'Gi = (Gi + (/- 2 ).( 5)52 and (2727 i)"(G2 - (5) = (G 2 - 5)9^ - 4727i[l] • <59. (6.11) 


Now using relations (16.61) and (j6.11l) the proof is straightforward by the induction. The inequalities n < p + 4 
for the function Giu and n < p + 3 for the function G 2 U are justified by the following observation. The 
right-hand sides of the expressions (I6.7l) - (j6.10p involve the functions ( 7271 )^( 1 ] for either j < n—2 (function 
Gilt) or j < n — 1 (function G 2 U). Note that we do not need to apply ( 7271 )" to 1 in the left-hand sides of 
(16.811 and (j6.10l) . According to Remark IT^ the function ( 727 i)'’[l] is well defined for all j < p + 2. Hence 
the condition n — 2 < p + 2, i.e., n < p -I- 4, is sufficient for the function Giu and the condition n — 1 < p -I- 2, 
i.e., n < p + 3, is sufhcient for the function G 2 U. □ 

Remark 6.6. In the above corollary, the condition q G C^[—b,b] is used only to ensure the existence of the 
functions (2727i)l[l], j < p + 2. The corollary holds whenever the corresponding number of the functions 
(27271 ) 1 [ 1 ] are well-defined independently on the smoothness of q. 

Now we show that the result of Corollary 16. 51 can be inverted in the following sense. Supposing that for a 
function r]{x) the functions ( 7271 )'^ [ 7 ], j = 1,2,... are well defined, we are going to show that the preimages 
of the function rj under the operators Gi and G 2 possess corresponding numbers of derivatives. 

Let u G C^[—b,b] be a given function and denote 7 := Giu. Applying Gf^ to the first relation in (| 6 . 6 I) 
we obtain 

2 G-S 27 i^ = Gf\G 2 u' - u'(0)) =u'- u'(0)G-i[l] 

= u' - u'(0) = dGf^r] - dGf^r]{0). 

It follows from the definition of the operator Gi that 77 ( 0 ) = u(0) and r]'{Q) = u\0) hu{0). Hence 

9 Gj"^ 77 ( 0 ) = m'(0 ) = r]'{0) — hr]{0) and we obtain the first inverse commutation relation 

G^^ [ 2727117 ] =9GC^77-r7'(0)-b/i?7(0). ( 6 . 12 ) 

In order to obtain the second commutation relation consider u G Gq [—b, 6 ] and denote rj = G 2 U. Applying 
Gf^ to the second relation in (j 6 . 6 l) we obtain 

2Gf^j2liV = Gf^{Giu' + f ■ u'(0) - 2u'(0)) = u'+ u'(0) - 2u'{Q)Gf^[\], (6.13) 
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where we used the inverse of the equality Gi[l] = /. It follows from the definition of the operator G 2 that 
77 '(0) = u'{0), hence (16.131) can be written as 

2Gr'7272^ = dG^^r,+ 7j'{0) - 277'(0)Gr'[1], 


or 

Gr' [272717 + 27 '( 0 )] = dG^^r, + 77'( 0 ). ( 6 . 14 ) 

Since the operator Gi is a bijection on G^[—b,b], the commutation relation (j6.12l) holds on the whole 
G^[—b,b]. Moreover, it means that for any 77 £ G^[—b,b] (which is equivalent to the condition that there 
exists 72717 and it is a continuous function), the preimage Gi^r] is a continuously differentiable function. 
Similarly, the commutation relation ( 16.1411 holds on Gq [—b, b] and allows one to conclude that the preimage 
G^^77 is also a continuously differentiable function. Since G2[l] = 1 , the commutation relation ( 16.1411 can be 
extended onto G^[—b,b] in the following way 

Gi^ [27271(7 - 7(0)) + 27 '( 0 )] = 5 G ”^77 + 77 '( 0 ). ( 6 . 15 ) 

Suppose now that for 77 £ Gg[—6, 5 ] there exists (27271)^77 and is a continuous function. Denote v := 
2727177. We obtain from ( 16 . 121 ) and ( 16 . 141 ) that Gi^v is continuously differentiable and 

G^\ = Gr'( 2727 i 7 ) = dG^^r, + rj'iO) - 2r^'{0) ■ Gr'[l]. 

If additionally Gj"^[l] is continuously differentiable (or equivalently, by Corollary 16.41 there exists 27271 [ 1 ]), 
then G^^77 is twice continuously differentiable and 

G^-'(( 27271 )^ 7 ) = 277'(0) • 5GrMl] - (2727 i7)'(0), (6.16) 

where we have used that 72717(0) = 0 . The equality ( 16 . 121 ) for 77 = 1 reads as 

G ^-1 (27271 [ 1 ]) =5Gr'[l] + / 7 , 

hence we can write (| 6 . 16 l) as 

G2-'((2727i)\) = d^G^^rj - 2 r^'{ 0 ) (G^"'(27271W) - h) - (2727177)'(0), 

or 


G2^ ((2727i)"7 + 27(0) • 2727i[1]) = d^G^^ + 2hrj\0) - (2727i7)'(0). (6.17) 

Similarly we obtain from (16.141) that Gi^rj is twice continuously differentiable function and that 

GC^ ((27271)^7 + 2(2727i7)'(0)) = d'^G^^r] + (2727i7)'(0). (6.18) 

In order to obtain the higher order differentiability of the preimages consider the following functions 

Mo = 7 ^^0 = 7-7(0) (6.19) 

U2n+1 = 27271 [W2n], l’2n+l = ‘^l2ll[v2n] + 2 m 2 „( 0 ), (6.20) 

M2n+2 = 27271 [7i2n+l] + 27i2„+i(0), M2n+2 = 27271 [7;2n+l], 71 = 0,1,... (6.21) 


We suppose that the recursive process stops when we cannot evaluate the operator 7271 anymore, i.e., the 
obtained function is not continuously differentiable. Note that the function Un is a linear combination of the 
functions (27271 )'^ 7 , j < n and the functions (27271)-^ [1], j < n —2. The function Vn is a linear combination 
of the functions ( 2727 i)'^ 7 , j < n and the functions ( 27271 )^[1], j < n (however if 7 ( 0 ) = 0 the function 
(27271 )"[ 1 ] does not participate). 

Starting from the relations (16.121) . (|6.14|) . (16.171) and (I6.18|) by induction we obtain the following result. 

Proposition 6.7. Let a function 7 £ C^[—b,b] be given and the functions Uj and Vj be well defined by the 
formulas (j6.19l) " (|6.21|) forj < n. Then the preimages Gj”^7 and Gf^rj are n times continuously differentiable 
functions and the following relations hold, 

d-Gf\ = Gf^[un] - m'„_ 2 ( 0 ), d-Gf^v = G^-'K] - 2 /i<_ 2 ( 0 ) + <_i( 0 ), 
d^Gf^p = - m;_i( 0 ), 5"Gr^ = G^^K] + m'„_i( 0 ) - /im„_i( 0 ). 
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if n is even, 
if n is odd. 


( 6 . 22 ) 

(6.23) 
































The existence of the functions Uj and Vj can be guaranteed by imposing smoothness conditions on q and 
rj, c.f., Remark 14.21 and Corollarv l6.5l The following corollary holds. 

Corollary 6.8. Suppose that q € C^[—b^h] and let a function rj € C'"[—6,6] be given. Then € 

G"*[—6,6], where m = min{n,p + 4}. 

Let a function Tj £ CQ[—b,b] be given. Then Gf^r] £ C^[—b,b], w;6ere m = min{n,p + 3}. 

The transmutation operators Gi and G 2 together with Corollary 16.81 allow us to prove a result similar 
to Proposition 15.31 under different assumptions on the potential q. Suppose that in some neighborhood of 
X = 0 the potential q possesses continuous derivatives of all orders up to the order N. Then both functions 
gi{x) = ^ + jJq q{s) ds and 52 ( 3 ^) = \ Jq 9(s) ds participating in Theorem l2.17l are A^+1 times continuously 
differentiable at a; = 0 and there exist generalized derivatives 71 ( 7271 )^ 51 . 2 , 5 < iV in a neighborhood of 
zero. By Theorem 14. 101 the following representations hold 

T + 7 / q{s)ds = +o(a;^+^) and i f q{s) ds = +o(a;^+^). (6.24) 

2 4 Vo ^0 4 Jo n\ 

By Corollary 16.81 the functions gi := Gf^gi and 52 := Gf^g 2 are TV + 1 times continuously differentiable at 
cc = 0. Since the operators Gf^ and Gf^ are bounded and map the functions c„ and s„ into 2"“^a;" (see 
Proposition 16.IL the following Taylor formulas are valid 


pn-l„ .^ra 

5 i(a;) = oo + V -^— + o(a;^+^) 


n—1 


nl 


f[+} on-lu „rt 

and g2ix) = V -^— + o(a;^+^). 


n—1 


nl 


Due to the definition of the operator G the functions gi and 52 coincide (up to some additive constants) 
with the functions (j){x) + ip{x) and (fip'J ~ J’i^) from (15.11) . respectively. Repeating the proof of Proposition 
15.31 and using the uniqueness of the Taylor coefficients one can obtain the following result. 


Proposition 6.9. Let q £ G^[—c, c], where [— c, c] C [—6,6] and let the coefficients a„, 
defined by (I6.24|) . Then 


r ^ 2 n — 6271, 


n = 1,.. 


+ 1 
2 


bn, n < N + 1 be 


7 Convergence rate estimates for the analytic approximation of 
transmutation operators 

In this section we establish relations between the smoothness of the potential q and the decrease rate of 
£ 1,2 from Theorem 12.171 as functions of N. As a result, we establish convergence rate estimates for the 
approximations Kn of the integral kernel K. 

Theorem 7.1. Let q £ G^[— 6, 6]. Then there exists a constant Cp and a sequence of the approximate kernels 
Kn of the form (I2.28|) such that 


ma.x_\K{x, t) - KN{x,t)\ < (7.1) 

Proof. Consider the functions gi{x) = ^ j Jq q{s) ds and 52 ( 2 ;) = j fj q(s) ds and their preimages gi := 
Gf^gi and 52 := Gf^g 2 . Since gi G G^+^J—6,6] and 52 G Gq'*'^ [—6, 6], by Corollary 16.81 we have that 
5 i G 6 , 6] and 52 G Gq'*’^]— 6 , 6]. By Jackson’s theorem [51 Chap.4, Sec.6] there exist a constant Cp 

and sequences of polynomials {pN}Neti and {qN}Neti such that 

niax |gi(a;) -pn{x)\ < and inax 152 ( 2 ;) - qN{x)\ < 

x^[-b,by ' l\p+^ x^[-b,by ' 

and additionally polynomials qN satisfy 5 Ar( 0 ) = 0. 
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Applying the operators Gi and G 2 we obtain that for each N there exist such complex numbers 
,..., ,..., that the constants Si and S 2 in (j2.26l) and (12.271) satisfy 


£1 < 


CpIlGill 

jyp+i 


and £2 < 


Cp\\G2\\ 

iVP+1 ’ 


hence 


c.f., m Theorem 5.1]. 


max_|K(a;,t) — KN{x,t) \ < c(ei +£ 2 ) < 


ATP+l ’ 


□ 


Remark 7.2. Let a function g S G^[—b,b] be given. Then its sequence of polynomials of the best uniform 
approximation {prtjnGN satisfies 


max \g(x) — pn{x)\ =0 


1 


N —>■ 00 , 


_/VP+i J ’ 

see, e.g., flSl Thm. VIII]. Henee one can obtain a slightly stronger result than (17.11) . namely that 

1 


m.a-K_\K.{x,t) — KN{x,t) \ =0 


_/VP+i f ' 


(7.2) 


For the estimate dUD to hold does not need to be a continuous function, it is sufficient that the 
derivative be bounded or that q^P~^^ G Lip[—6,5]. We do not enter into further details in the present paper. 

The order of the convergence rate in Theorem l7.1l is close to optimal as the following partial inverse result 
shows. Before we need the following lemma. 

Lemma 7.3. Suppose that for a function u there exist all generalized derivatives ( 2727 i)-^m for j < k on 
a segment [—d,d\ C [—5,5], the last generalized derivative {2^27i)^u € G^[—d,d\ and q € G'^[—d,d], where 
k, £ and m are some non-negative integers. Then the function u possesses minj/c + £, m + 2} continuous 
derivatives on the segment [—d, d]. 

Proof. The general solution of the equation 27271 jy = u can be written in the form 

,(x) = c,/W + «/W^ J^ + /(X)^ 

Hence ii u G G^ and / € C" then g € The result of the lemma follows by induction taking into 

account that for q S d, d] we have / G d, d]. □ 


Theorem 7.4. Suppose that there exist a constant p gN, constants d > 0 and Cp > 0, and for each N there 
exist such complex constants ..., b[^\ ..., 5^^ that the constants £1 and £2 in (12.261) and (12.271) 

satisfy 

max{£i,£2}< (7-3) 

Then q G GP{—b, 5). 


Proof. Again, consider the functions gi{x) = ^ + z fo 92 {x) = 7 q{s) ds and their preimages 

gi := Gf^gi and g 2 ■= Gf^g 2 . Due to Proposition 16.II the preimages 


\n^0 / \n^l 

are polynomials of degree less or equal to N. Since the operators Gf^ and Gf^ are bounded, the inequality 
(m implies that the functions gi and <72 can be approximated by sequences of polynomials in such a way 
that the approximation errors decrease at least as 0{N~p~^~^). Due to the inverse approximation theorem 
(see, e.g., [El Theorem 31]), 51 G GP+^{-b,b) n C'-^P+^'>G][-b,b] and 32 G CP+^(-5, 5) D 5]. 
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Let [— (i, d] C {—b,b) and assume that q € C™‘[—d,d] but q ^ d,d] where m < p. In such case 

applying Corollary 16.51 to the function gi on the segment [—d, d] we obtain that there exist generalized 
derivatives ( 27271)-^51 of all orders up to A: = min{p + 1, to + 4}. Due to the assumption to < p we have that 
k > m + 2. Lemma 17.31 provides that gi £ d, d], contradicting the assumption q ^ [—d, d] as 

q = 9'i- □ 

Remark 7.5. There is another possibility to obtain convergence rate estimates similar to those of Theorems 
and \7.4\ It is based on the following simple observation: for each fixed x, the expression 

2[Ar/2] + l 2[(Ar+l)/2] 

KN{x,t) := aoUo{x,t) + ^ anU 2 n-l{x,t) + ^ bnU 2 n{x,t), (7.4) 

n—1 n—1 

where [•] denotes the integer part function, is a polynomial in t whose degree is less or equal to N. All 
generalized wave polynomials Un which are not included in (17a contain terms with powers of t higher 
than t^, and expression (Tfa contains exactly 2 terms with the power t^, these terms are (up to nonzero 
multiplicative constants) (po{x)t^ and (pi{x)t^. The x-derivative of (17.41) possesses the same properties, the 
terms with t^ are (pQ{x)t^ and ip({x)t^. 

The mentioned properties are sufficient to verify that the Cauchy problem 

= (7,5) 

K{b,t) = pi{t), Ki{b,t) =p2{f), 

where pi and p 2 are some polynomials of degree not greater than N, possesses a unique solution in the region 
0 < |t| < a: < 6 and this solution has the form (17.41) . 

Now the analogue of Theorem \7.1\ can be obtained if one approximates K(b, t) and Ki{b, t) by polynomials 
and utilizes continuous dependence (see, e.g., ]S2f ) of the Cauchy problem (Il3 on the initial data. The 
analogue of Theorem [13 can be obtained by utilizing continuous dependence of the solution of the Goursat 
problem (12.121) . (12.131) on the data on the characteristics and reasoning as in the proof of J21[ Theorem 3.5]. 

The described method does not require any inverse operators however the order of the convergence rate 
estimate which can be obtained is less by 1 than those provided by Theorem \71\ The reason is that for 
q £ 6] one has K(b,f) £ CP'^^[—b,b] while Ki{b,t) £ CP[—b,b] and hence the convergence rate of the 

polynomial approximations of Ki has estimates which are one order less than those for the function K. 


8 Transmutation operators restricted to the half segment 

For the solution of spectral problems it may be convenient to make use of the known initial values (12.211) 
and (12.221) . For that, one can consider equation (12.41) restricted to the segment [0,&]. The solutions c{uj,x) 
and s{uj,x) can be calculated once one knows the integral kernel K.{x,t) on the domain i?i = {(x, t) : 0 < 
\t\ < X < b}, and it can be deduced from the integral equation (I2.16P that the integral kernel K is uniquely 
determined on Ri by the parameter h and values of the potential q on [0, b] and does not depend on the 
values of qioi x < 0. That is, let the potential q be given on [0, b]. One can consider any sufficiently smooth 
continuation of q onto the segment [—b, b], the resulting integral kernel of the transmutation operator on Ri 
does not depend on the continuation. However the analytic approximation of the transmutation operator 
given by Theorem 12.171 does, as well as the inverse operators and the preimage of the integral 

kernel k. 

Note that the functions pfe, tpk, c^, s^, k £ No, are completely determined on [0,6] by the values on 
[0, b] of a particular solution / and instead of considering the approximation problems (I2.26P and (12.271) 
on the whole segment \—b,b] one can solve them on the half-segment [0,6]. It is clear that the coefficients 
ao,..., On, bi,... ,bN determined as solutions of the approximation problems (12.261) and (12.271) on the half¬ 
segment [0, 6] would not necessary solve them on the whole segment [—6, 6], it depends on the choice of the 
continuation of the potential q. We already verified in various numerical examples (see [27]) that the method 
based on restrictions of the approximation problems (12.261) and (|2.27l) to the half-segment works. In this 
section we briefly justify the main results for this restricted case. 

First we prove the main approximation theorem. 
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Theorem 8.1. Let the complex numbers uq, ■ ■ ■, on and bi,... ,bN be such that the inequalities (12.261) and 
(12.271) hold for any x S [0,&]. Then the kernel K(a;,i) is approximated by the function (12.281) in such a way 
that for any (x, t) G Ri the inequality holds 

|K(x,<) - K]s[{x,t)\ < C{£i +£ 2 ), 

where 

C = 3/o(6\/M), M := max|g(x)|, 

[0,b] 

and Iq is the modified Bessel function of the first kind. 

Proof. Consider the difference = K — K]\f and the corresponding function H^{u, v) = K^(u + v,u — , 

Q < u + V < b. It is the unique solution of the Goursat problem (12.101) , (12.111) and due to the inequalities 
(I2.26P and (I2.27P the functions i7®(u,0) and 77®(0,u) are bounded by ei +£ 2 . The Goursat problem (12.101) . 
(12.111) is equivalent to the following integral equation 

pU pV 

H%u,v) = H^{u,0) + H^{0,v)-H^0,0)+ / q{u'+ v')H{u',v') du'dv', 

Jo Jo 

whose solution can be estimated by the successive iterations method, see, e.g., m Proposition 18]. Note 
that |i7®(0, 0)1 < £1 + £ 2 , hence the successive iterations provide that 

|i7®(u, u)| < 3 (£i + £ 2 ) - = 3(£i + £ 2 )Io{\/M\uv\) < 3(£i + £ 2 )Io{b'/M) . 


Approximations of the derivatives c'{ui,x) and s'{uj,x) require considering additionally a transmutation 
operator for the Darboux-associated equation 

— u" + qD(x)u = LO^U, (8.1) 


where 


to = 2 ( y ) - q. 


The function 1// is a particular solution of (18.11) for a; = 0 satisfying (l//)(0) = 1 and (l//)'(0) = —h. 
Consider the transmutation operator Ti/f for (18.11) . where the subindex 1// means that the operator Ti/f 
is such that Ti/j[l] = Ijf. Denote the integral kernel of Ti/f by K^. We refer the reader to ^7}, [12] for 
further details. 

As was shown in m, the pair of operators T/ and Ti/f allows one to write the derivatives of the 
solutions of (I2.15P in the following form 

f 

c'(uj,x) = —c{uj,x) — ujTi/ f[suiujx], 

f 

s'(a;, x) = + Ti/j^[cosa;a;]. 

Recall that the generalized wave polynomials for the operator Ti/f have the form 


Vo=tpo{x), V2m-l{x,t) 


m ^ < 

{^]'Pfn-k{x)t'", 

even A;=0 '' ' 


V27n{x,t) 


m y V 

(^)'(pm-kix)t’'- 

odd fc=l J 


Denote by c„ and s„ the corresponding traces, c.f., (I2.24|) and (|2.25l) . 

It was proved in [I?] Theorem 6.4] that if one constructs an approximation of K via Theorem 12.171 then 
the expression 

N N 


Kd.n = -boVo{x,t) - ^ bnV 2 n-l{x,t) - ^ a„U 2 „(x,t), (8.2) 

n—1 n—1 

where bo := ao, can be used as an approximation to the integral kernel Kb. We present a simpler proof of 
this result, also suitable for the case of the half-segment. 


25 






























Theorem 8.2. Let the complex numbers ag, ■ ■ ■, on and bi,... ,bN be such that the inequalities (I2.26P and 
(12.271) hold for any x S [0, b]. Then the kernel 'Kd(x, t) is approximated by the function (18.21) in such a way 
that for any (x, t) € Ri the inequality holds 


where 


JCni^jt) — < C{eiMi + (ei + £2){2MiM2b + Mi + 1)^, 


C = 3Io{bVM), 


M := max |g_D(x)|, 

[0,b] 


Ml := max 

[0,b] 


1 

fix) 


M 2 := max \f'{x)\ 

[o,&] 


(8.3) 


and Iq is the modified Bessel function of the first kind. 


Proof. Consider the functions 


h 1 ^ 

9^ix) = -2~^7 QDis)ds+ '^bnCnix), 
•^0 n=0 

92 ( 2 ;) = -/ qpjs) ds + a„s„(x). 

^ Jo _ 1 


(8.4) 


(8.5) 


Note that it is sufficient to show that the functions gi and 52 are small enough to derive (18.31) from Theorem 
18.11 applied to equation (18.11) . Consider also the functions 


h 1 ^ 

5i(a;)=2 + 4/ qis)ds-'^GnCnix), 

•^0 n=0 

1 /.a; N 

92 ix) = - g(s)ds-V5„s„(x), 

^ Jo _ 1 


the left-hand sides of the inequalities (12.261) and (12.271) . 
One can easily verify that the following equalities hold 


[%D{s)ds) =-fdU^ + \ [%is)ds), 


f ^ j,^m — ^(Sm—1 “1“ — — j,dfSrnj 

fd — S>fji = ?7i(c77^_i + = —OfCrri') 

f9jCo = J^f^o = 0, 

where for convenience we denoted Sg = Sq = 0. Hence 


m = 1,... 
m = 1,... 


and 


f^j9i = -jdf92 


fdj92 = -jdf9i- 


Inverting the equalities ()8.6I) and (j8.7l) one obtains that 


5i(a:) = 


Cl 


1 


fix) fix) Jo 


fis) 


Cl . s , 52(0) 

- 92ix) + + 


92is) 

fis) 

2 


fix) 


fix) fix) Jo 


ds 

fis)g2is)ds 


( 8 . 6 ) 

(8.7) 


( 8 . 8 ) 
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and 


f'is)gi{s)ds, 


(8.9) 


52(2;) 


C2 

fix) 


9i{x) + 


gi(0) 

fix) 



where ci and C 2 are some constants. Substituting x = 0 into ira . (EH) and (18.41) . (j8.5l) one can find that 
Cl = bo — ^ = Go — ^ = —gi(0) and C 2 = 0. Finally, using the inequalities |( 7 i(x)| < ei and \g 2 ix)\ < 62 , 
X £ [0,6], one obtains the following estimates 


|gi(a:)| < eiMi + e 2 ( 2 MiM 2 |a:| + Ml + 1) and |g 2 ( 2 :)| < ei( 2 MiM 2 |a;| + Mi + 1), 
which combined with Theorem 18.11 finish the proof. 


□ 
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